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Abstract: 

Qh! We show how a difference in the correlation length longitudinally and transversely, with 
respect to the jet axis in e + e~ annihilation, arises naturally in a model for Bose-Einstein 
correlations based on the Lund string model. In genuine three-particle correlations the 
difference is even more apparent and they provide therefore a good probe for the longitudinal 
stretching of the string field. The correlation length between pion pairs is found to be rather 
independent of the pion multiplicity and the kaon content of the final state. 
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1 Introduction 



The Hanbury-Brown-Twiss (HBT) effect (popularly known as the Bose-Einstein effect) cor- 
responds to an enhancement in the two identical boson correlation function when the two 
particles have similar energy-momenta. A well-known formula [p]] to relate the two-particle 
correlation function (in four momenta pj,j = 1,2 with q = p\ — p?) to the space-time density, 
p, of (chaotic) emission sources is 

where TZ is the normalised Fourier transform of the source density 

/ pjx) expjiqx) 

n{q) = JpWTx (2) 

The commonly used event generators HERWIG and JETSET are based upon classical 
stochastical processes and do not include HBT-effects (although Sjostrand, in JETSET, 
has introduced an ingenious method to simulate any given distribution by means of a kind 
of mean-field potential attraction between the bosons in the final state). 

In this letter we will further investigate some features of the methods developed in |3j (an 
extension of to multi-boson final states). We will show that the model predicts, due to 
the properties of string fragmentation, a difference between the correlation length along the 
string and transverse to it. In practice this means that if we introduce the longitudinal and 
transverse components of the vector q (defined with respect to the thrust direction) then 
we obtain a noticeable difference in the correlation distributions. This becomes even more 
noticeable when we go to the three-particle HBT effect (which was predicted in ||) because 
in this case even more of the longitudinal stretching of the string field becomes obvious. 
Finally we will investigate the influence of the kaon and baryon content of the states on the 
HBT effects between the pions. 



2 Longitudinal and transverse correlation lengths 

The starting point of our Bose-Einstein model [3] is an interpretation of the (non- 
normalised) Lund string area fragmentation probability for an n-particle state (cf Fig P 

n 

dP( Pl ,p 2 , . . . ,p n ) = \[Ndp 3 5(p) - m*)6(£ Pj ~ Ptot) exp(-bA) (3) 

i 

in accordance with a quantum mechanical transition probability containing the final state 
phase space multiplied by the square of a matrix element Ai. In and in more detail in [[| 
a possible matrix element is suggested in accordance with (Schwinger) tunneling and the 
(Wilson) loop operators necessary to ensure gauge invariance. The matrix element is 

M = exp(m - b/2)A (4) 

where the area A is interpreted in coordinate space, k is the string constant (phenomenolog- 
ically k ~ 1 GeV/fm) and b ~ 0.3 GeV/ fm is the decay constant. Note that the parameter 
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Figure 1: The decay of a Lund Model string spanning the space-time area A. The particles 1 
and 2 are identical bosons and the particle(s) produced in between them is denoted by I. The 
two possible ways, (...,1,1,2,...) and (. . . , 2, 1, 1, . . .), to produce the state are shown and 
the area difference between the two cases, AA, is shaded. The two neighbouring vertices of 
the state with the two identical bosons are denoted by a and h, and the transverse momenta 
of the quarks produced in the neighbouring vertices are ±k_L a and ±k^ b , respectively. 

b is much smaller than k. From now on we will, as is usual in the Lund model, go over to the 
energy momentum space. Then the area A — > 2k 2 A, while b — > b/2n 2 , as explained in 

The transverse momentum properties are in the Lund model taken into account by means 
of a Gaussian tunneling process. In this way the produced (gg)-pair in each vertex will 
obtain ±k_i_ and the hadron stemming from the combination of a q from one vertex and a q 
from the adjacent vertex obtains = kj_j+i — k±j. 

In case there are two or more identical bosons the matrix element should be symmetrised 
and in general we obtain the symmetrised production amplitude 

M = Y,M V (5) 

v 

where the sum goes over all possible permutations of the identical particles. The squared 
amplitude occurring in Eq (|3|) will then be 

JETSET will provide the outer sum in Eq (||) by the generation of many events but it is 
evident that the model predicts a quantum mechanical interference weight, w-p, for each 
given final state characterised by the permutation V: 

1 _ 2Re(M-pM* P ,) 
wv = l+Y, , M 2 J ' (7) 
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In the Lund Model we note in particular for the case exhibited in Fig [I], with two identical 
bosons denoted 1 and 2 having a state I in between, that the decay area is different if the 
two identical particles are exchanged. It is evident that the interference between the two 
permutation matrices will contain the area difference, AA, and the resulting general weight 
formula will be 

AA 



w v = i + 



cos 



2 




cosh ( + A(Ek 

where A stands for the difference between the configurations described by the permutations 
V and V' and the sum is taken over all the vertices. In our MC implementation of the weight 
we replace the string constant k in the transverse momentum generation with the default (in 
JETSET) transverse width, 2cr 2 (which is of the order of k). The calculation of the weight 
function for n identical bosons contains n\ — 1 terms and it is therefore from a computational 
point of view of exponential- type. We have in |3j introduced approximate methods reducing 
it to power-type instead and we refer for details to this work. 

We have seen that the transverse and longitudinal components of the particles momenta 
stem from different generation mechanisms. This is clearly manifested in the weight in 
Eq (|8|) where they give different contributions. In the following we will therefore in some 
detail analyse the impact of this difference on the transverse and longitudinal correlation 
lengths, as implemented in the model. 

In order to understand the properties of the weight in Eq (§) we again consider the 
simple case in Fig [I]. The area difference of the two configurations depends upon the energy 
momentum vectors pi,p2 and pi and can in a dimensionless and useful way be written as 

A A 

— = 5 P 5x L (9) 
2k 

where 5p = P2 — pi and 5xl = (St; 0, 0, Sz) is a reasonable estimate of the space-time 
difference, along the surface area, between the production points of the two identical bosons. 

In order to preserve the transverse momenta of the particles in the state (1,7,2) it is 
necessary to change the generated kj_ at the two internal vertices around the state / during 
the permutation, i.e. to change the Gaussian weights. Also in this case we may write a 
formula similar to Eq (|9|) for the transverse momentum change: 



A(£k 



2 ^ 
-Li' 



2k 



Sp ± Sx ± (10) 



where 5p ± is the difference p_i_ 2 — P_li and <5xj_ = (k^ b — (— kj_ a ))/K. The two neighbouring 
vertices of the state (1,1,2) ((2,7,1)) are denoted by a and b and k^ b + kj_ a corresponds 
to the states transverse momentum exchange to the outside. Therefore <5x_|_ constitutes a 
possible estimate of the transverse distance between the production points of the pair. 

For the general case when the permutation V' is more than a two-particle exchange there 
are formulas similar to Eqs @ and ([T0| ) although they are more complex (and the expressions 
do not vanish when only two of the exchanged particles have the same energy momentum). 
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It is evident from the considerations leading to Eqs (^j) and (|T0| ) that only particles with a 
finite longitudinal distance and small relative energy momenta will give significant contribu- 
tions to the weights. We also note that we are in this way describing longitudinal correlation 
lengths along the colour fields, inside which a given flavour combination is compensated. 
The corresponding transverse correlation length describes the tunneling (and in this model 
it provides a damping chaoticity). 

The weight distribution we obtain is discussed in (and with varying kaon and baryon 
content also below). It is strongly centered around unity although there are noticeable tails 
to both larger and smaller (even negative) weights. The total production probability is, 
however, positive and we find negligible changes in the JETSET default observables (besides 
the correlation functions) by this extension of the Lund model. 



3 Results 

Two-dimensional Bose-Einstein correlations in e + e~ annihilation have been analysed at lower 
energies than LEP by the TASSO collaboration H|. Although they find that their data is 
compatible with a spherically symmetric correlation function they conclude that at least one 
order of magnitude of more data is required to obtain more detailed information. With the 
large statistics available from LEP we have therefore generated gg-events at the Z° pole to 
investigate the properties of our model. Short-lived resonances like the p and K* are allowed 
to decay before the BE-symmetrisation, while more long-lived ones are not affected. 

We have analysed two-particle correlations in the Longitudinal Centre-of-Mass System 
(LCMS). For each pair of particles the LCMS is the system in which the sum of the two 
particles momentum components along the jet axis is zero, which of course also means that 
the sum of their momenta is perpendicular to the jet axis. The transverse and longitudinal 
momentum differences are then defined in the LCMS as 



Qt = V (P*2 - Pxl) 2 + (p y 2 - Pyl) 2 (H) 
<?L = \Pz2 ~Pz2\ 

where the jet axis is along the z-axis. 

We have taken the ratio of the two-particle probability density of pions, P2, with and 
without BE weights applied as the two-particle correlation function, R 2 

D / \ p2w(Pl,P2) 

R2{PuP2) = - f r- 12 

P2{P1,P2) 

and the resulting function is shown in Fig ^| It is clearly seen that it is not symmetric in 
and q± and in particular that the correlation length, as measured by the inverse of the width 
of the correlation function, is longer in the longitudinal than in the transverse direction. 
This difference remains for reasonable changes of the width in the transverse momentum 
generation. For comparison we have also analysed events where the Bose-Einstein effect has 
been simulated by the LUBOEI algorithm implemented in JETSET In LUBOEI the BE 
effect is simulated as a mean-field potential between identical bosons which is spherically 
symmetric in Q. Analysing only the initial particles and particles stemming from short-lived 
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Figure 2: The ratio Rz{qLiq±) of the number of charged pion pairs having relative four- 
momentum components qi and q± with and without Bose-Einstein weights applied. The 
sample consists of particles which are either initially produced or stemming from short-lived 
resonances. 

decays results for the LUBOEI events in a correlation function with identical transverse and 
longitudinal correlation lengths. The correlation lengths are in agreement with the source 
radii input to LUBOEI. Using all the final pion pairs, after all decays, in the analysis results 
in a small decrease in the transverse correlation length and of course a large decrease in 
the height for qi ~ q±_ ~ 0, while the longitudinal correlation length is rather unaffected. 
The pions from long lived decays affect the correlation lengths in the same way both for our 
model and for LUBOEI. 

In [||] it is shown that our model gives rise to genuine three-particle correlations. We will 
in this letter continue to investigate three-particle correlations and we will in particular use 
our knowledge of the different contributions to the weight function to study the genuine 
higher order correlations. We will also exhibit how the genuine higher order terms in the 
weight function mainly clusters particles in the longitudinal direction. 

The total three-particle correlation function is in analogy with Eq (O) 



T>"( \ P3w{Pl,P2,P3) nQ , 

Pa\Pi,P2,P3) 

To get the genuine three-particle correlation function, R%, the consequences of having two- 
particle correlations in the model have to be subtracted from R 3 . To this aim we have 
calculated the weight taking into account only configurations where pairs are exchanged, 
w . In this way the three-particle correlations which only are a consequence of lower order 
correlations can be defined as 

R 3 {puP2,P3) = — t r 14 

Pz\Pl,P2,Pz) 
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Figure 3: R 3 {Q) and R 3 {Q) are shown in the left figure, while the figure to the right shows 
R 3 (Q). The existence of genuine three-particle correlations is apparent. 

The genuine three-particle correlation function, R 3 , is then given by 

R 3 = Rl - R' 3 + 1 (15) 

We have analysed R3 in one dimension as a function of the kinematical variable 

Q = ^Q 2 12 + Q 2 13 + Q 2 23 with Ql = -( Pl - Pj y (16) 

and in two dimensions we have used the following variables calculated in the L CMS for each 
triplet of identical bosons 

QL = \l qll2 + QL13 + QL23 Wlth llij = (Pzi-Pzj) 2 (17) 

q± = v^iiT+^iiT+^i^ with ilij = (p±< - p^) 2 

where the ^-axis is along the jet axis. In Fig ^| the correlation functions R 3 (Q), R 3 (Q) and 
i?3(Q) are shown, and the existence of genuine three-particle correlations in the model is 
clearly exhibited. 

This way of getting the genuine correlations is not possible in an experimental situation, 
where one has to find other ways to get a R 3 reference sample. We have suggested one possible 
option in || and the results in this letter are in agreement with the conclusions of that 
investigation. In the present analysis the contribution to the correlations from higher order 
configurations in the weight calculation is apparent. We note that R 3 flattens out earlier, 
i.e. for lower Q- values than R 3 . This means that the genuine three-particle correlations 
have a longer correlation length compared to the consequences of lower order correlations. 
Performing the same analysis in two dimensions in the LCMS for each triplet results in the 
R^QliQa.) distribution shown in Fig |j. The effect of the higher order terms is to pull the 
triplets closer in the longitudinal direction while the transverse direction is rather unaffected. 
This suggests that higher order correlations are more sensitive to the longitudinal stretching 
of the string field. 

We have also studied the correlation length for pion pairs as a function of the final charged 
multiplicity and the kaon content of the state. Within statistical errors which are relatively 
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Figure 4: The ratio Rz{qL-,<I±) of the number of triplets of charged pions with and without 
Bose-Einstein weights applied. 

large we see no dependence on either the charged multiplicity or the number of kaons. Since 
one might suspect that events with many pions are premiered by the re-weighting the average 
baryon and kaon content of the events have been investigated. We find that the changes 
of the average multiplicity of different kaon species as well as of the average multiplicity of 
protons and neutrons in the final state are much smaller than the experimental errors as 
summarised in 01. 
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